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Stability  of  Rectilinear  Geostrophic  Vortices 
in  Stationary  Equilibrium 

by 

Louis  Bauer  and  George  K.  Morikawa 

Courant  Institute  of  Mathematical  Sciences 
New  York  University,  New  York,  N.  Y.  10012 

Abstract 

A hemispherical  barotropic  model  of  global-scale 
atmospheric  motions  is  studied  in  a planar  approximation 
using  three  Kq  (Bessel)-vortices  of  equal  strength  y and 
length  scale  /c*1  to  represent  the  sub-tropical  semi- 
permanent high  pressure  systems.  In  equilibrium  these 
anticyclonic  vortices  are  equally  spaced  on  a latii ade 
circle  of  radius  a.  One  KQ-vortex  of  strength  -yo  is  used 
to  simulate  the  polar  (cyclonic)  vortex;  is  chosen  so 
that  the  equilibrium  configuration  remains  stationery 
with  respect  to  the  rotating  earth.  Over  the  range  of 
the  length  parameter  0 < (jca)  _<  5>  two  distinct  ste  ble 
regimes  are  found:  1)  a " soft"-vorcex  regime, 

0 _<  (< a)  < (/c^r*a)  = 1.44763516  and  2)  a "hard" -voi  tex 
(climatological)  regime,  k > <cr-  The  nonlinear  pert  of 
this  stability  investigation  has  been  carried  out  mainly 
by  numerical  computation. 
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I.  Introduction 

Rectilinear  geostrophic  (Bessel)  vortices  were  first  used  by 
Stewart  [1,2],  to  study  the  periodic  properties  of  the  motion  of 
the  semi-permanent  atmospheric  pressure  systems.  He  found  that 
the  period  of  these  oscillations  is  of  the  order  of  magnitude  of 
years.  He  represented  the  subtropical  anticyclones  by  sin  equilib- 
rium configuration  of  three  vortices  of  equal  strength  y and 
length  scale  /c-1  spaced  equidistantHy  on  a latitude  circle  of 
radius  a;  but  this  equilibrium  has  a slow  westerly  drifting  motion. 
Another  equilibrium  configuration,  which  he  considered  with  an 
additional  cyclonic  vortex  fixed  at  the  pole  to  counteract  the 
drift,  was  found  to  be  unstable  [2]. 

In  this  paper  we  study  the  motion  of  geostrophic  vortices 
about  an  equilibrium  configuration  in  which  the  polar  (or  center) 
vortex  is  free  to  interact  with  the  three  anticyclonic  vortices. 

The  strength  yQ  of  the  polar  vortex  is  chosen  so  that  the  equilib- 
rium configuration  is  stationary  with  respect  to  the  rotating 
earth.  The  center-vortex  length  parameter  <c  = < is  usually  taken; 
but  the  effect  of  different  / ic  is  discussed  briefly.  We  esti- 
mate that  the  climatological  range  of  the  dimensionless  length 
parameter  to  be  2.5  < (<a)  < 5»  The  inclusion  of  a free  polar 
vortex,  albeit  a rather  simplified  representation  of  the  totality 
of  vorticity  concentrations  at  the  higher  latitudes,  introduces 
two  additional  degrees  of  freedom,  resulting  in  a physically 
meaningful  model  of  atmospheric  motions  on  a hemispherical  scale. 
The  present  study  is  an  outgrowth  of  the  recent  work  of  Moi Jkawa 
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and  Swenson  [3]  where  the  possibility  of  finding  stable  motions 
for  this  stationary  equilibrium  configuration  is  indicated.  For 
the  nonlinear  problem  we  say  that  the  motion  is  unstable  if  the 
center  vortex  goes  outside  of  the  equilibrium  circle.  The  compre- 
hensive stability  investigation  carried  out  in  Ref.  [3]  shows  that 
the  usual  normal-mode  type  of  linear  stability  analysis  is  not 
sufficient  to  determine  stable,  oscillatory  motions,  especially 
due  to  the  stationarity  condition.  The  stationarity  condition 
introduces  two  additional  zero  eigenvalues  which  yield  algebraic 
instabilities  in  addition  to  two  pairs  of  eigenfrequencies  ±0-, 

?.rd  tOg*  This  implies  that  the  full  nonlinear  equations  must  be 
studied  in  order  to  reveal  and  ascertain  the  stable  motions. 
Numerical  computations  are  made  for  this  purpose. 

The  geostrophic  vortex  representation  is  a particular  aspect 
of  the  barotropic  model  [4]  of  atmospheric  motion.  The  principal 
approximations  are:  l)  heating  and  dissipation  are  neglected, 

2)  hydrostatic  balance  replaces  the  vertical  momentum  equation, 

J)  geostrophic  balance  replaces  the  horizontal  momentum  equations 
and  4)  the  Coriolis  parameter  is  constant,  fixed  at  a particular 
latitude,  for  example,  at  the  pole.  Our  results  indicate  that, 
for  the  very  large-scale  pressure  systems  of  planetary  dimensions, 
free  long-period  oscillatory  motions  are  possible. 

In  addition,  a more  comprehensive  numerical  study  has  been 
made  over  the  full  range  of  the  length  parameter,  0 < (<a)  _<  5* 

A remarkable  phenomenon  is  found  with  respect  to  nonlinear  stabil- 
ity over  this  ra.nge  of  (< a).  Normalizing  by  setting  a = 1,  the 
stable  ranges  of  k are  divided  into  two  parts:  1)  a low  k.  range. 


0 _<  k < 1.4476  and  2)  a high  k range,  k > 1.4476  (including  the 
climatological  range).  The  approximate  value  * = 1.4476  corre- 
sponds precisely  to  the  frequency  ratio  = 2,  at  which  value 

the  motion  is  unstable  for  any  Initial  perturbation.  Both  stable 
ranges  are  of  physical  interest.  A comparison  of  the  behavior  of 
the  stable  nonlinear  interacting  motions  in  the  two  stable  regimes 
shows  that  the  low  k vortices  act  like  "soft"  vortices  and  the 
high  k vortices  act  like  relatively  "hard"  vortices.  The  low  k 
vortices,  especially  for  k « 1,  are  exceptionally  stable  for  a 
circle-vortex  initial  displacement;  whereas  the  hard  vortices  in 
the  upper  climatological  regime,  e.g.  * = 5,  are  very  stable  for 
center-vortex  initial  displacements. 


II.  Equations  of  Motion  and  Condition  for  Stationary  Equilibrium 

We  study  the  interacting  motion  of  geostrophic  vortices 
where,  initially,  the  vortices  are  positioned  near  a stationary 
equilibrium  configuration  which  consists  of  N (anticyclonic ) 
vortices  of  equal  strength  7 equally  spaced  on  a (latitude)  circle 
of  radius  a,  and  one  (cyclonic)  vortex  of  strength  yQ  at  the 
center  (north  or  south  pole)  of  the  circle;  the  configuration 
(with  N =■  3 ) which  we  study  is  shown  in  Fig.  1.  A single  vortex 
is  given  by  the  streamfunction 

(1)  ip  = yKQ(<r)  , r = (x2+y2)1/? 

where  positive  7 corresponds  to  anticyclonic  (clockwise)  motion 
and  Kq  is  the  modified  Bessel  function  of  the  second  kind  (zeroth 
order).  The  velocity  components  are 

(2) 

which  also  correspond  to  the  conditions  for  geostrophic  motion. 

The  2(N+1)  kinematic  equations  of  motion  [3]  for  interacting  geo- 
strophic vortices  in  Cartesian  coordinates  are  (for  k = 1,  2,  ...,N) 
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(5) 

(6) 


K, (p  . ) 

lvpoi ' 

Poi 

Kl^oj) 

Poi 


where  the  dot  means  time  derivative;  is  the  Bessel  function  of 
the  second  kind  (first  order);  and 

Pki  = *hxk-x1)£  + (yk  -yj)2]172 
pko = *t(Xk  -x0)2  + (yk  -y0)2]1/£! 

Poi  - ^b„^1)!t(y0-yi)s]1/! 

The  following  dimensional  normalizations  have  been  made:  Lengths 

are  normalized  by  the  circle  radius  a;  time  is  normalized  by 
(a  /y);  and  the  center  vortex  strength  yQ  is  normalized  by  7.  in 
most  of  the  work  described  in  this  paper,  one  length  scale  K-1  is 
associated  with  ach  of  the  vortices.  However  we  have  carried  out 
in  addition  some  calculations  with  another  length  scale  for 
the  center  vortex  in  order  to  simulate  different  values  of  the 
C^iiolis  parameter  at  the  pole  than  at  the  latitude  circle;  for 
these  computations,  (3)  to  (6)  are  modified  by  replacing  k.  by  * 
in  those  terms  having  the  strength  factor  7q  of  the  center  vortex. 

The  general  equilibrium  condition  [3]  is  uniform  rotation  of 
the  circle  vortices  with  angular  velocity 

(?a)  0=0  +y  a) 

0 1 o o 

with 


I 


5 


o. 


-hi 


N 


and 


7 4 r—  akiKl^aki^ 


ojq  = tcK^ic) 


where 


aki  = " cos  Ww  )1 


1/2 


ki 


w.  . = w.  - w . = w, 
ki  k l k 


-w°=f  (K-i) 


N 


The  superscript  denotes  the  initial  equilibrium  position.  From 
(7a)  the  condition  for  stationary  equilibrium  is  n = 0,  or 


(7b) 


T'o  * -°o/a>o 


For  the  nonlinear  numerical  computations  we  use  one  of  the 
integrals  [5]  of  the  motion  [Ref.  3,  Appendix  B,  Eqn.  ( B9 ) ] : 


(8) 


N 


k=0 


Vk 


r.  = const, 


(8)  is  used  to  monitor  the  error  accumulation, 


r 


III.  Linearized  Stability  Analysis 


For  N _>  3 the  2(N+1)  linearized  equations  of  motion  (constant 
coefficients)  for  the  perturbed  positions  ) 

from  equilibrium  are 

( 9 ) rf 1 ) + y A,  . r ! 1 ) + 5 B,  . ( 0 [ 1 ) - df1^) 

k iWTT^i  1 iWTT^l  kl  1 k J 

/~(l)  . 0 ~(l)  O \ 

+ Vo11;  sin  wk'yo  cos  V j 0 


N N 

(10)  9^1)  +(Q  + C0+70D0)41)+  } Ckiri1)+^  \±d± 

k o o o k i^j=1  ki  i i^JZ1  Tti  i 


(1) 


" 7nDjx1^  cos  w£+%1}  sin  w°)  = 0 


(ID  ^1)+R0T^^1)  sin  wJ-cd0  ^ e{1}  cos  w°  - (0  + Eo)y(1)  =0 
(12)  ^-D.^rp)  cos  sinwJ  + fO  + E^-O 


Where 


cki 


°u 


= (-<  /2 ) sin  w.  .K  (a,  . ) 

' ki  o'  ki 

= (<2/2)(  (1 +cos  wki)K0(cki)  + [2K1(cki)/cki]) 

= U‘/2)( (i  . cos  wki)K0(oki)  + [ 2Ki( aki )/ak±^ ) 

N 

- wki)Ko(0kl>  -2  C0E  wki[Kl(oki)/oki11 

= <[2(1 -cos  wki)]V?  , Wlr.  = w°,  = (2tt/N)  (k-i ) 


’ki  "ki 


0) 


= *%(<)  , Dq  = <2(K0(<)  + [K1(<  )/<]  1 , Eo  = (N/2)<2Kj(<) 
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The  position  of  a circle  vortex  is  described  in  polar 

.<vf  i ) l ) . 

coordinates  and  the  position  of  the  center  vortex  (x^  ,y^  ')  is 


described  in  a Cartesian  frame  of  reference  rotating  with  the 
equilibrium  angular  velocity  O.  For  0=0,  (x^^y^^ ) = (xj^^yj^  ) 


o 7 ^ o 


in  the  rest  frame.  For  no  center  vortex  -y  = 0 (O  ^ 0),  (9)  and 
(10)  reduce  to  equations  equivalent  to  those  studied  by  Stewart 
[1].  Since  (9)  and  (12)  have  constant  coefficients,  the  usual 


exponential  stability  analysis  is  possible  as  described  in  Ref.  3* 
In  the  remainder  of  this  section  we  summarize  the  results  of  these 
calculations  over  the  range  of  physical  parameters  of  interest. 

First  we  review,  correct  and  extend  the  calculations  for 
Stewart's  model  [1]  with  N = 3 and  -y  = 0.  With  2N  degrees  of 
freedom,  there  are  six  eigenfrequencies  ±0^  of  which  two  are  zero 
eigenfrequencies  [5].  The  non-zero  eigenfrequencies  Oc  and  Oq1, 
where  Oo  is  the  angular  velocity  of  the  equilibrium  around  the 
circle,  are  shown  in  Table  I for  k = 3>  3»5  and  4.  Based  on 
Stewart's  estimates  this  range  of  c (or  (*ca))  corresponds  nominally 
to  a latitude  range  between  35°  and  20°.  We  note  that,  for  N = 3 
and  a 20°  latitude  circle,  the  equidistance  between  anticyclones 
is  about  10,000  km;  this  is  also  nominally  the  distance  between 
pole  and  equator  where  the  geostrophic  concept  tends  to  become 
inapplicable . 

Table  I.  Non-Dimensional  Eigenfrequencies  for  (N=3i  yc  =0) 


K 

°o 

°ol 

0 

£ 

1 — 1 

0 

c 

3 

0.01690 

0.05711 

3.38 

3.5 

0.00761 

0.02899 

3. 81 

4 

0.00340 

0.01440 

4.24 
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The  anticyclonic^  vortex  strength  7 is  estimated  from  the  tangential 
velocity  relation  using  (1)  and  (2) 

(13)  |V( r ) | = 7<K1(<r) 

Following  Stewart's  estimate,  the  length  scale  k-1  = 2000  km.  is 
assumed;  for  r = 2000  km.,  the  characteristic  velocity  is  estima- 
ted to  be  of  the  order  of  10  meters  sec-1;  then  the  strength 
6 2 -1 

7 = 2.87  x 10  km  day  . Now  the  oscillatory  periods  t^  corre- 
sponding to  the  non-dimensional  frequencies  are  obtained  by 

(14)  ti  = 27ra2/~yOi  = 8.756(<a)2/oi  days  cycle-1 

where  a is  the  distance  between  pole  and  latitude  circle.  The 
oscillatory  periods  t and  t ^ corresponding  to  Oq  and  o ^ in 
Table  I are  given  in  Table  II. 

Table  II.  Oscillatory  Periods  for  (N=5;  7q=0) 


U a) 

t 

0 

*01 

days  (years) 

days  (years) 

3 

4665  (12.8) 

1580  (5-78) 

3.5 

4 

41200  (115) 

9750  (26.7) 

For  5 _<  (/ca)  _<  4,  the  time  t for  an  anticyclone  to  drift  once 
around  the  latitude  circle  ranges  between  15  and  113  years. 

Now  we  introduce  a free  cyclonic  polar  (center)  vortex  in 
order  to  nullify  this  drifting  motion.  First  we  investigate  the 
exponential  stability  for  all  N when  the  stationarity  condition 
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O = 0 is  imposed.  From  Ref.  3 we  find  that  a circle  vortex  is 
exponentially  unstable  for  (N  _>  4;  o =••  0)  and  any  value  of  < > 0. 
This  is  illustrated  in  Fig.  2 where  the  dashed  curve  for  (N  = 4; 

O = 0)  is  on  the  stippled  unstable  side  of  the  neutral  stability 
curve  for  N = 4.  For  (N  = 2;  O = 0)  the  center  vortex  is 
exponentially  unstable  as  shown  in  Fig.  2.  Only  for  (N  = 3;  O = 0 
is  the  stationary  equilibrium  curve  on  the  exponentially  stable 
side  of  the  corresponding  neutral  stability  curve.  With  2(N+l) 
degrees  of  freedom,  there  are  eight  eigenfrequencies  ±0^  of  which 
four  are  zero  eigenfrequencies  [6],  The  non-zero  eigenfrequencies 
Oj  and  and  the  value  of  yQ  from  the  stationarity  condition  (18) 

are  shown  in  Table  III  for  < = 2.5,  3,  3»5,  4,  4.5,  and  5*  Here 
the  same  length  scale  k ^ has  been  associated  with  both  center  and 
circle  vortices. 

Table  III.  Non-Dimensional  Eigenfrequencies  for  (N=3,  0 = 0) 


^0 

°1 

°2 

°2/°1 

°0/°l 

-0.2010 

0.1498 

0.6619 

4.417 

0.2478 

-0.1403 

0.07694 

0.5057 

6.573 

0.2197 

-0.09780 

0.03837 

0.3771 

9.829 

0.1984 

n n^firiQ 

n m Q7 0 

1 1 1 70 

VJ • U-LO  ( c. 

-LH  • [ c. 

Xo  J-O 

-0.04737 

0.008988 

0.1978 

22.01 

0.1679 

-0.03294 

0. 004260 

0.1399 

32.85 

0.1564 

The  oscillatory  periods  t-^  and  t ^ corresponding  to  the  non- 
dimensional  eigenfrequencies  O-^  and  o2  in  Table  III  are  again  cal- 
culated using  (l4)  in  the  range  5 < (< a)  _<  4 for  = 2000  km. 


* 


But  in  the  extended  range  2.5  ^ (<a)  _<  3,  we  choose  k~  - 2500  km 
and  t = 10.945 ( k a ) /oi;  and  in  the  extended  range  4 _<  (>ca)  < 5, 
we  choose  <_1  = 1500  km  and  ti  = 6.567 ( «:a)2/oi . Then,  for  the 
entire  range  2.5  _<  (xa)  _<  5>  the  range  of  the  corresponding  lati- 
tude circles  lies  between  18°  and  56°.  These  oscillatory  periods 
are  given  in  Table  IV. 


Table  IV.  Oscillatory  Periods  for  (N=3;  0*0) 


t„ 


Ua) 


Lat . , 
Deg. 


days  (years)  [ days  (years) 


kT1  =2500  km 


456  (1.25) 


1280  (3.51) 


103  (0.283) 


195  (0.534) 


1 = 2000  km 


3 

36° 

1024  (2.81) 

156  (0.427) 

3.5 

27° 

2795  (7.66) 

284  (0.779) 

4 

18° 

7484  (20.5) 

508  (1.39) 

kT1  =1500  km 


36° 

5613  (15.4) 

381  (1.04) 

29° 

14790  (4o.5) 

672  (1.84) 

22° 

38540  (106) 

1170  (3.21) 

Table  IV  shows  that,  over  the  length  parameter  range  2 .5  (ica)  < 5, 
the  long  period  range  is  I.25  yrs  _<  t1  < 106  yrs  and  the  short 
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period  range  is  0.283  yrs  < tg  < 3.21  yrs.  The  effect  of  changing 
the  vortex  length  scale  by  500  km.  is  indicated  for  two  values  of 
(<a)  = 3 and  4.  For  these  dynamically  balanced  vortex  motions, 
shorter  periods  are  associated  with  vortices  having  larger  length 
scales  (smaller  <). 

Since  the  four  zero  eigenvalues  yield  algebraic  instabilitie 
proportional  to  t or  V , the  validity  of  these  linearized  calcula- 
tions for  (N  = 3;  O = 0)  depends  on  the  study  of  the  nonlinear 
equations  (3)  to  (6)  by  numerical  computation. 


IV.  Numerical  Computation  of  Nonlinear  Motions 


A.  Type  of  Initial  Perturbations 

We  follow  the  procedure  established  in  Ref.  3 arid  restrict 
the  initial  perturbation  displacement  away  from  equilibrium  to 
only  one  vortex,  either  the  center  vortex  or  one  of  the  circle 
vortices.  Specifically,  four  different  types  of  initial  perturba- 
tions were  considered:  (a)  Displacement  of  the  center  vortex, 

xo  and  (b)  displacement  of  the  3-rd  circle  vortex, 

^(O)  and/or  y^  ^(0)  (see  Fig.  l).  In  this  section  the  super- 
script ( 1 ; notation  is  used  to  indicate  all  perturbation  displace- 
ments, not  necessarily  in  the  linear  range,  but  also  in  the  stable 
nonlinear  range  (as  in  Ref.  3 the  perturbation  parameter  e has 
been  included  in  the  displacement  to  simplify  the  notation).  From 
these  single  initial  perturbations  plus  symmetry  considerations, 
we  are  able  to  construct  the  range  of  stable  initial  displacements 
for  center  and  circle  vortices.  In  the  future  it  would  be  of 
interest  to  consider  combined  initial  displacements  of  more  than 
one  vortex,  as  well  as  random  initial  perturbations  for  one  run  or 
an  ensemble  of  runs. 

B.  Stationarity  Condition  for  Finite  Initial  Perturbations 

The  condition  for  stationary  equilibrium  is  given  by  (8). 
However  for  various  finite  initial  perturbations  especially  for 
large  initial  displacements,  even  though  the  motion  is  stable  in 
the  exponential  sense,  there  will  be  an  appreciable  drifting  motion 
unless  70  is  changed.  We  have  found  that  this  drift  is  essentially 
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a uniform  rotating  motion  and  can  be  eliminated  by  changing  the 
strength  ratio  y / /.  Thus  this  drifting  motion,  arising  from 
large  initial  perturbations,  may  be  interpreted  as  an  algebraic 
instability  which  is  removable  by  correcting  7/7.  Therefore  a 
corrected  7Q(/ca^)  i^  computed  so  that  the  circle  vortices  remain 
approximately  stationary  in  the  mean,  where  a^  is  the  distance 
between  the  i-th  circle  vortex  and  the  center  vortex.  As  a first 
order  correction  we  compute  an  average  distance  a by 


(15) 
and 

(16) 


a = (1/5)  ^ ai 
= -O0Ua)/w0Ua) 


In  many  cases  is  sufficient  to  nullify  the  drifting  motion. 
Although  a complete  elimination  of  the  drifting  motion  is  possible 
for  all  stable  cases,  we  have  usually  not  felt  the  necessity  for 
making  finer  corrections.  In  fact  a small  drifting  motion  often 
allows  us  to  discern  the  form  of  the  circle  vortex  trajectories 
more  easily  when  they  do  not  overlap  completely  after  several 
periods  of  motion.  Of  course,  when  the  center  vortex  becomes 
unstable  and  migrates  outside  of  the  circle,  the  circle  vortices 
undergo  a drifting  motion  with  angular  velocity  of  the  order  of  Oq. 


C.  Numerical  Methods 

Originally  a Runge-Kutta  [7]  method  was  used  to  solve  the 
system  of  eight  [2(N+l);N  = 3]  ordinary  differential  equations  (5) 
to  (6)  describing  the  interacting  motion  of  geostrophic  vortices. 
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This  procedure  is  the  following:  A basic  time  interval  size  h is 

o 

chosen.  The  integration  is  performed  over  each  interval  with  a 
spacing  hQ  and  = hQ/2.  If  the  maximum  relative  difference  in 
the  two  integrations  is  not  less  than  e « 1,  we  halve  the  time 
step  width  again,  h 2 = h^/2,  etc.,  and  repeat  the  procedure  over 
the  subinterval.  When  two  consecutive  integrations  agree  within  e 
the  integration  is  carried  out  over  hQ  with  the  required  sub- 
interval size.  At  the  next  interval  hQ  the  procedure  is  repeated. 
We  observe  (and  print  out)  the  actual  number  of  subintervals 
required  at  each  interval.  If  the  number  of  subintervals  becomes 

excessive  (8  or  16),  we  reduce  the  basic  interval  size  h and  redo 

o 

the  problem.  We  experimented  with  values  of  e of  10_R,  10-7  and 

-6  -6 

10  , and  found  that  e = 10  was  satisfactory  for  integrations 

over  1000  intervals,  i.e.  the  change  in  the  final  solution  for 

-6  -8 

e = 10  and  10”  did  not  reach  the  fourth  significant  digit  in 
the  worst  case.  It  is  in  fact  doubtful  whether  a much  smaller 
value  of  e is  advantageous  because  a large  number  of  subintervals 
leads  to  accumulation  of  round-off  errors.  The  standard  total 
time  T for  any  k was  taken  to  be  2tt/60  where  6Q  is  defined  by  (16). 
Usually  the  1000  time  steps  was  in  the  range  v/ 2C)Q  < T < ir/C)  . 

Eventually  when  the  accuracy  of  the  computations  was  assured, 
we  used  a modified  Adams  method  [8]  for  our  integrations,  but 
without  re-subdividing  each  interval.  The  results  are  substanti- 
ally the  same  as  the  ones  obtained  by  the  Runge-Kutta  method  and 
appreciably  faster  (about  25  seconds  per  run  on  the  CDC  6600 
instead  of  1^5  seconds). 


15 


D.  Numerical  Results  for  the  Climatological  Regime 


1.  Range  of  Stable  Initial  Displacements 

In  the  range  of  the  length  parameter,  2.5  < k < 5>  stable 
nonlinear  motions  are  found  over  a wide  range  of  initial  perturba- 


tion displacements  of  center  or  circle  vortices.  For  k < 2 
(02/01  < 3),  the  magnitude  of  stable  initial  perturbations  de- 
creases rapidly  and  approaches  zero  as  < -*•  1.4476  — 2) . 

This  is  illustrated  in  Fig.  5 for  initial  radial  displacements  of 
a circle  vortex  and  in  Fig.  4 for  initial  displacements  of  the 
center  vortex.  These  figures  show  that  a second  stable  regime 
exists  in  the  small  length-parameter  range  0 < k < 1.4476 
(1  < n < 2);  this  complementary,  soft-vortex  regime  will  be  de- 
scribed later  in  Section  IV-E.  In  the  climatological  hard-vortex 
regime  comparison  of  Figs.  3 and  4 shows  that  the  stability  bound- 
ary for  center-vortex  initial  displacements  is  smooth  and  monotonic 
with  respect  tu  the  length  parameter  k (Fig.  4);  however,  the 
stability  boundary  for  circle-vortex  initial  displacements  shows 
some  non-monotonic  wave-like  structures  — particularly  evident  in 
the  range  < ic  < 3»75.>  which  is  in  the  middle  of  the  estimated 
climatological  regime.  The  rib-like  lines  radiating  out  of  the 
(C-axis  are  lines  of  constant  frequency  ratio,  n (=  nQ  = f^/C^  on 
the  /c-axis),  and  n increases  with  increasing  k from  the  value 

n = 1 at  k = 0.  A more  detailed  view  of  the  above  wave-like 
o 

structures  is  found  by  xooking  at  the  stability  boundary  in  th^ 
two-dimensional  configuration  space  of  center-  and  circle-vortex 
initial  displacements  for  several  values  of  Figure  5 shows 
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both  center-  end  circle-vortex  stability  boundaries  for  < = 3,  3*5 
and  4;  Figs.  6 and  7 show  the  circle-vortex  stability  boundaries 
for  k = 3.2  and  3.5,  respectively;  and  Fig.  8 shows  the  center- 
vortex  stability  boundary  for  k.  = 4.  As  shown  in  Figs.  5 and  8, 
all  the  center-vortex  stability  boundaries  are  slightly  cusped 
equilateral  triangles  with  smooth  faces  directed  towards  the  three 
circle  vortices.  The  vertices  of  the  center-vortex  stability  tri- 
angle are  greater  than  0.46  from  the  center  for  < = 5 (Fig.  4)  and 
still  increasing  for  k > 5;  from  geometrical  considerations,  an 
apparent  limit  is  the  vertex-center  distance  < 0.5*  In  the  clima- 
tological (hard-vortex)  regime  the  center-vortex  stability  region 
is  substantially  larger  than  the  circle-vortex  stability  region. 

In  the  soft-vortex  regime  0 < k < 1.4476,  this  inequality  is 
reversed.  Figure  6 for  k = 3.2  shows  the  most  irregular  structure; 
in  fact  there  are  ragged  peninsulas  and  stability  islands;  the 
stability  boundary  is  only  approximately  symmetrical  with  respect 
to  the  radial  line  through  the  equilibrium  position  — we  see  no 
reason  to  expect  exact  symmetry.  Contour  lines  of  constant  n are 
shown;  the  line,  n = 5,  goes  through  the  stability  islands..  In 
our  rather  comprehensive  numerical  calculations,  we  may  still  have 
missed  a few  islands.  For  < = 3*5  (Fig.  7)  there  is  only  one 
slender  island,  paralleling  the  main  stability  boundary  and  lying 
approximately  symmetrical  with  respect  to  the  radial  line  through 
the  equilibrium  position;  the  contour  line,  n = 5#  passes  through 
the  island  and  the  stability  boundary  is  almost  coincident  with 
the  contour  line,  n = 6.  Figure  3 and  comparison  between  Figs.  6 
and  7 show  that  the  major  wave-like  break-di  ^n  of  the  stability 
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boundary  for  large  amplitude  perturbations  of  a circle-vortex 
occurs  between  the  frequency  ratio,  5 < n < 6.  Even  without  addi- 
tional numerical  computations,  the  structure  of  these  irregular 
stability  boundaries  in  the  three-dimensional  parameter  space 
[/c,*x^1^(0),y^1^(0)]  can  be  visualized  for  3 < < < 3.75.  For 
example,  as  the  length  parameter  is  changed  from  K =3.2  to 
k = 3»5 > the  small  islands  in  configuration  space  disappear  and 
the  peninsulas  shrink  into  one  slender  island;  for  < > 3.65,  this 
island  disappears.  In  additior  to  this  major  irregularity  in  the 
stability  boundary  for  5 < n < 6,  we  find  smaller  stability- 
boundary-break-down regions  for  8 < n < 9 and  11  < n < 12,  etc. 

't  is  remarkable  that  this  rather  symmetrical  configuration,  with 
only  eight  degrees  of  freedom  (four  Bessel  vortices)  and  with 
vortices  having  smooth,  monotonic  stream  functions  and  velocity 
distributions,  yields  such  a variety  of  structure  in  its  stability 
characteristics  for  single- vortex  initial  displacements. 

2.  Oscillatory  Periods  and  Trajectories  of  Nonlinear  Motions 

The  linearized  stability  analysis  described  in  Section  III 
shows  that  for  (N  = 3>0  = 0)  and  k.  > 0 there  are  two  pairs  of  non- 
zero eigenfrequencies,  o1  and  02  and  corresponding  periods  (cf. 
Tables  III  and  IV).  We  now  interpret  these  periods  in  terms  of 
the  vortex  trajectories  for  small  initial  displacements  equal  to 
0.001  with  k = 2:  In  Fig.  9a,  x^1^(0)  = 0.001;  in  Fig.  9b, 

x-^(O)  = -0.001;  in  Fig.  9c,  y^(0)  = 0.001;  and  in  Fig.  9d, 

xo1^(°)  ~ 0.001.  These  trajectories  are  enlarged  500  times.  The 
vortex  trajectories  for  negative  initial  displacements  are 
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obtainable  from  the  vortex  trajectories  for  positive  initial  dis- 
placements by  rotating  each  vortex  trajectory  about  its  own  equi- 
librium position  by  tt  (or  one-half  revolution,  cf.  Figs.  9a  and 
9b).  In  each  case,  the  center-vortex  trajectory  has  counter- 
clockwise loops  progressing  clockwise  about  the  origin  with  a 
single  excursion  taking  place  after  four  loops;  and  subsequently 
recurrence  takes  place.  In  all  cases  the  normalized  recurrence 
time  (or  period)  Oq/Or  = 0.2864,  which  agrees  with  the  linearized 
long  period  00/01  (Oq  = 0.0803)  for  * = 2.  In  the  time  tR  = 2rr/0  , 
each  circle  vortex  makes  a counterclockwise  excursion  around  a 
closed  trajectory  near  its  equilibrium  position.  In  the  linearized 
analysis  the  frequency  ratio  nQ  = {0^0^  = (t-j/t^;  for  K = 2, 
nQ  = 3.000.  We  define  the  short  period  as  tR/n;  from  the  center- 
vortex  trajectory  for  < = 2 with  small  initial  displacement,  n is 
equal  to  (four  counterclockwise  loops)  minus  (on'  clockwise 
excursion).  Generally  n is  an  integer  only  for  certain  discrete 
values  of  <,  where  n increases  for  increasing  K.  Between  two  con- 
secutive integer  values  of  n,  say  n,  and  n0,  there  will  be  two- 
excursion  recurrence  for  K corresponding  to  a frequency  ratio 
(nl  + 1/2),  three-excursion  recurrence  for  k's  corresponding  to 
f re  }uen^y  ratios  (n^  + 1/3)  and  (n1  + 2/3),  etc.;  for  irrational 
values  of  n there  will  be  no  recurrence  at  all.  From  the  center- 
vortex  trajectory  with  recurrence,  the  general  algorithm  for 
rational  values  of  n is: 

n = ((nujn^er  °f  loops)  - (number  of  excursions  in  one  recurrence)] 

(number  of  excursions  in  one  recurrence) 
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This  algorithm  for  determining  the  frequency  ratio  n is 


applicable  for  all  initial  displacements  with  recurrent  trajec- 
tories in  the  entire  range  of  stable  motions.  In  the  climato- 
logical range  (n  > 2)  for  fixed  k. , n decreases  monotonically  from 
its  linear  value  as  the  initial  displacement  is  increased  to  the 
stability  boundary  (see  Figs.  3 and  4).  The  dependence  of  the 

A 

normalized  recurrence  time  O /Od  on  n for  center-  and  circle- 
vortex  initial  displacements  is  shown  in  Fig.  10  for  k = 2.5,  3, 

4,  4.5  and  5*  For  each  <,  the  upper  curve  is  for  negative 
initial  radial  displacements  of  a circle  vortex  x^^(o)  <0, 
y^(0)  = 0 and  the  lower  curve  is  for  negative  initial  displace- 
ments of  the  center  vortex  x^(0)  < 0,  y^(o)  = 0.  All  the 
curves  are  apparently  linear;  and  specifically  for  x'  ^(0)  < 0, 
the  short  period  tR/n  = const.  = tg  * t-j/n  with  less  than  5% 
deviation.  For  each  k,  the  curves  for  all  other  types  of  initial 
displacements  are  bracketed  by  these  two  curves.  The  dependence 
of  the  vortex-strength  ratio  |yj_1  on  n for  center-  and  circle- 
vortex  initial  displacements  is  shown  in  Fig.  11  for  k = 2.5,  3, 
3«5>  4 and  5.  For  each  >c,  the  upper  curve  is  for  positive  initial 
radial  displacements  of  a circle  vortex  xi1^(0j  > 0 and  the  lower 

curve  is  for  negative  initial  radial  displacements  xi^(O)  < 0. 

3 

The  curves  for  all  other  initial  displacements  are  bracketed  by 
these  two  curves,  including  the  curve  for  negative  initial  dis- 
placements of  the  center  vortex  x^  ^(O)  < 0 which  is  shown 

/N 

explicitly.  The  normalized  recurrence  times  t„/t,  and  the  vortex- 
strength  ratios  | -y Q | - 1 for  integer  values  of  n near  the  stability 
boundary  are  given  in  Table  V for  circle-vortex  initial  displace- 
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Eli] 


[x^1)(o),y^1}(o)] 


Vh 


(-0.126, o) 
( 0.220,0) 
( 0.240,0) 
(-0.158  0) 
(-0.160, 0) 
( 0.215,0) 
(-0.160,  0) 
(-0.165,0) 
( 0.216,0) 


4.62 

5.68 

8.47 

8.95 

12.58 

18.07 

17.73 

24.86 

39.35 


0.783 

0.795 

0.763 

0.699 

0.686 

0.789 

0.706 

0.657 

0.796 


Table  VI.  Recurrence  Times  for  Center-Vortex 
Initial  Displacements  x^^O) 


vfc 


K. 

[41)(0),y^1)(o)] 

2.5 

(-0.28',  0) 

5.3^ 

3 

(-0.360, 0) 

7.74 

3.5 

(-0.400,  0) 

11.54 

4 

4.5 

(-0.410,0) 
(-0. 44^,0) 

17.00 

25^ 

5 

(-0.450,0) 

33.29 

R'  "1 


0.566^ 

0.396 

0.315 

0.276 

4 

0. 182 ' 
0.148 


j Estimates  from  Figs.  4,  10  and  11. 
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ments  x-^ 1 ^ ( 0 ) and  in  Table  VI  for  center-vortex  initial  displace- 
ments x^^(O).  The  real  recurrence  times  tR  can  be  calculated 
from  the  linearized  values  of  the  long  period  t^  given  in  Table  IV. 
For  example,  for  [/c,  n/x^1  ^ ( 0), y^1  ^ ( 0)]  = ( 5>  5; -0.450,  0), 
tR  = 0. 15  x 106  yrs  = 16  yrs  and  the  short  period  tR/n=3«2  yrs  = t,  . 

The  center-vortex  trajectory  for  small  initial  perturbations 
is  confined  in  a circular  annulus,  which  changes  into  a triangular 
annulus  for  large  initial  displacements  approaching  the  stability 
boundary.  For  both  center-  and  circle-vortex  perturbations,  the 
center-vortex-trajectory  loops  become  cusped  for  large-amplitude 
motions  near  instability.  When  the  motion  becomes  unstable,  the 
center  vortex  goes  outside  of  the  circle.  Several  typical  trajec- 
toiies  are  shown  in  the  ra  ige  2.5  < < < 5-  Fig.  12  shows  the  tra- 


jectories for  a small  initial  tangential  displacement  of  a circle 
vortex,  [<,n,-x^1^(0),y^1^(0)]  = ( 2.66, 5; 0,  0. 001 );  comparing  to 
Fig.  9,  the  center-vortex  trajectory  is  seen  to  loop  out  for  a 
y^^O)  displacement  and  loop  in  for  an  displacement.  Two 

series  of  recurrent  trajectories  near  the  stability  boundary, 
including  unstable  trajectories  outside  the  stability  boundary,  are 
shown  in  Figs.  13  and  l4  for  < = 3*2:  1)  In  Fig.  13a, 

[(c,n;x|^(0),yj^(0)]  = (3-2, 6;  0.196,0)  and  = 0.860;  in  Fig. 

13b  with  [x-^^Oby^1^  (0)]  = (0.215,0),  the  center  vortex  becomes 
unstable  after  approximate?. v four-seventh  of  the  recurrence  time 
for  n = 6 (cf.  Fig.  13a);  in  Fig.  13c  [k,  n;  x^1  ^ ( 0),  y^ 1 ^ ( 0)] 

= (3.2,5/0.230,0)  and  tR/t1  = 0.789  and  the  initial  displacement 
lies  on  the  stability  island  on  the  x-axis  in  Fig.  6;  in  Fig.  13d 
with  [x^^Ohy^  (0)]  = (0.240,0),  the  center  vortex  becomes 
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unstable  after  approximately  a full  recurrence  time  for  n = 5 (cf. 

Fig.  13c).  And  2)  in  Fig.  l4a,  C«,n;x^1)(0),y^1)(0)] 

= (3 • 2, 5; -0. 120, 0.208 ) with  t^/t^  = 0.810  wherein  the  circle- 

vortex  initial  displacement  has  a large  tangential  component  and 

the  initial  position  lies  on  the  upper  peninsula  in  Fig.  6;  in 

Fig.  l4b,  [ic,njx^1^(0),y;^1^(0)]  = (3*2,9/2;-C. 130,0. 240)  and 
/\ 

t^/ti  = 0.787  wherein  the  initial  position  lies  near  the  stability 
boundary  on  the  upper  peninsula;  in  Fig.  14c  with  [x^1  ^ (Oby^1^  (0)] 
= (-0. 120,  0.235  )>  the  center  vortex  becomes  unstable  after  more 
than  two  recurrence  times  for  n = 9/2  (cf.  Fig.  l4b).  Comparison 
of  series  1)  (Fig.  13)  with  series  2)  (Fig.  l4)  shows  some 
features  in  common  as  well  as  contrasting  or  distinguishing 
characteristics  of  trajectories  for  radial  vs.  predominantly  tan- 
gential initial  displacements.  The  common  features  are:  the 

circle  vortices  make  the  same  number  of  excursions  about  their 
equilibrium  positions  as  the  center  vortex  about  the  configuration 
center;  usually  the  stable  trajectories  are  confined  within  the 
stability  boundaries  with  the  exception  of  initial  displacements 
on  islands  or  peninsulas;  near  the  stability  boundary  the  center 
vortex  tends  to  make  cusps  for  any  type  of  large  initial  displace- 
ment; there  is  bilateral  symmetry  of  the  stable  configurations 
(after  a recurrence  time)  with  respect  to  a radial  line  parallel 
or  nearly  parallel  to  the  x-axis;  the  unstable  center  vortex  first 
becomes  captured  by  one  of  the  three  circle  vortices,  but  then 
changes  partners,  progressing  in  a clockwise  order  about  the 
center.  In  the  least  unstable  of  the  three  unstable  cases  shown 
for  k = 3»2,  the  center  vortex  returns  to  the  center  of  the 
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configuration  for  a brief  time  (Fig.  l4c).  The  contrasting 
1 matures  are:  for  radial  initial  displacements,  the  center  vortex 

makes  inside  loops  about  the  center  and  the  circle  vortex  trajec- 
tories tend  to  loop  out  or  make  cusps  (Fig.  13);  while  for  pre- 
dominantly tangential  displacements  the  center  vortex  makes  out- 
side loops  (away  from  the  center)  and  cusps  while  the  circle  vortex 
trajectories  are  relatively  smooth.  Next,  an  example  of  recurrent 
trajectories  for  a center-vortex  initial  displacement  near  the 
stability  boundary  is  shown  in  Fig.  15  for  [ «c,  n; x£ 1 ^ ( 0),  y ^ ^ (0) ] 

= (2 4-J 0. 1955*  0)  and  t^/t^  =0.508.  The  center-vortex  trajectory 
is  confined  in  a narrow  annulus  so  that  a second-order  drift 
correction  of  -11.1$  is  required  in  order  to  resolve  the  frequency 
ratio  n = 4;  also  the  circle  vortex  trajectories  are  difficult  to 
resolve.  Finally  in  Fig.  16  for  [^n/xj^O^yj^O)] 

“ (4 «5* 15; 0. 2138, 0)  and  t^/t^  = 0.785,  an  example  is  shown  of 
recurrent  trajectories  for  a circle-vortex  initial  displacement 
near  the  stability  boundary  with  a relatively  large  K and  n. 

E.  Soft  Vortices  Compared  to  Hard  Vortices 
1.  Mainly  Stable  Motions 

The  behavior  of  the  soft  vortices  in  the  small-length- 
parameter  range  0 < K < 1.4476  (1  < n < 2)  is  compared  with  the 
characteristics  of  the  hard  (climatological)  vortices.  The  varia- 
tion of  the  stability  boundary  with  respect  to  k for  initial 
radial  displacements  along  the  x-axis  of  a circle  vortex  and 
center  vortex  have  already  been  shown  in  Figs.  3 and  4 respec- 
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tively.  Comparing  the  inclination  of  the  n-contours  (constant 
frequency  ratio  02/01)  for  large  displacements  in  the  soft-vortex 
regime  with  the  slope  of  the  n-contours  in  the  hard-vortex  regime, 
both  sets  of  contours  are  directed  away  from  the  instability 
point,  kcj,  = 1.4476  (n  = 2)  for  increasing  values  of  the  initial 
displacements.  The  numerical  calculations  show  that  in  the  range 
1 < /c  < 2,  which  straddles  the  instability  point,  the  stability 
boundary  itself  approximates  a node-like  contour,  n = 2;  this  is 
most  evident  for  negative  initial  radial  displacement  of  a circle 
vortex  (or  a positive  initial  displacement  of  the  center  vortex) 
on  the  soft-vortex  side  of  <cr  and  for  positive  initial  radial 
displacement  of  a circle  vortex  (or  a negative  initial  displace- 
ment of  the  center  vortex)  on  the  hard-vortex  side  of  k.  ; see 

cr 

Fig.  17  for  [ /c,  n; xj 1 ^ ( 0 ) , y-j 1 ^ ( 0 ) ] = (1.25,2; -0. 118,  0)  with 
tR/ti  = O.6175  and  a second-order  drift  correction  (approximately 
-1%)  was  not  made;  when  a drift  correction  is  incorporated  into  the 
run  shown  in  Fig.  17,  this  run  is  clearly  shown  to  be  unstable, 
i.e.  any  run  with  n = 2 is  unstable.  Figure  3 shows  that  with  a 
circle-vortex  initial  displacement  for  small  k.  (or  n approximately 
equal  to  one)  the  soft-vortex  configuration  is  appreciably  more 
stable  than  the  hard  vortex  configuration  (with  n > 2).  In 
contrast  Fig.  4 shows  that,  with  a center-vortex  initial  displace- 
ment, the  soft-vortex  configuration  is  less  stable  than  the  hard- 
vortex  configuration.  In  Fig.  18  the  stability  boundaries  are 
shown  in  the  two-dimensional  configuration  space  and  center-  and 
circle-vortex  initial  displacements  for  k = 0.75;  the  vertices  of 
the  triangle-like  stability  boundaries  are  directed  toward  each 
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other.  In  contrast  the  "vertices"  of  the  stability  boundaries  for 
center-  and  circle-vortex  initial  displacements  in  the  hard-vortex 
regime  shown  in  Fig.  5 are  directed  away  from  each  other. 

The  dependence  of  the  normalized  recurrence  time  Q /o~  on 

o R 

the  frequency  ratio  n for  some  center-  and  circle-vortex  initial 
displacements  is  shown  in  Fig.  19  for  soft  vortices  k = 0,  0.25, 
0.50,  0.75>  1»0  and  I.25.  For  each  < the  upper  curve  is  for 
positive  initial  radial  displacements  of  a circle  vortex 
x^^^(o)  > 0,  y^ 1 ^ ( 0 ) = 0;  and  the  lower  curve,  which  first  dips 
below  the  linear  curve,  is  for  negative  initial  radial  displace- 
ments of  a circle  vortex  x^^(o)  < 0,  y^ 1 ^ ( 0 ) = 0.  The  curves  for 
ol±1  other  initial  displacements  lie  between  these  two  curves, 
including  the  single  curve  for  all  (±)  center-vortex  initial  dis- 
placements which  is  shown  explicitly  for  two  xr's,  k = 0.75  and  1.0. 
In  contrast  the  corresponding  hard-vortex  curves  all  lie  below  the 
linear  curve  (see  Fig.  10 ). 

The  dependence  of  the  vortex-strength- ratio  ly0l-1  on  the 
frequency  ratio  n for  circle-vortex  initial  displacements  is  shown 
in  Fig.  20  for  ic  = 0,  0.25,  0.50,  0.75  and  1.0.  For  each  ic,  the 
upper  curve  is  for  positive  initial  radial  displacement  of  a 
circle  vortex  x^1^(o)  > 0,  y^^O)  = 0;  and  for  < = 0.75  and  1.0 
the  lower  curve  is  for  negative  initial  radial  displacement  of  a 
circle  vortex  x^ 1 ^ ( 0 ) < 0,  y^ 1 ^ ( 0 ) = 0.  Curves  for  other  initial 
displacements  are  bracketed  by  these  curves  which  have  a 
small  spread  in  the  soft-vortex  regime.  In  contrast  the  corre- 
sponding hard-vortex  curves  all  lie  above  the  linear  curve  (see 
Fig.  11).  Note  that  the  coordinates  are  much  enlarged  in  Fig.  20 
as  compared  to  Fig.  11. 
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A number  of  recurrent  trajectories  of  soft-vortex  motion  are 


shown  next,  mainly  near  the  stability  boundary.  In  addition 
1)  for  < = 0.125,  a run  is  shown  without  a second-order  drift 
correction  to  compare  with  the  corresponding  compensated  recurrent 
trajectories  without  the  drifting  motion;  and  2)  for  k = 0 the 
small  initial  perturbation  runs  are  shown  as  well  as  two  unstable 
runs  beyond  the  stability  boundary.  A comprehensive  series  of 
recurrent  trajectories  for  initial  radial  displacements  of  a circle 
vortex  for  < = 0.75,  0.50,  0.25,  0.125  and  0 follows:  In  Fig.  21, 

U,n;  x^1,(0),yj1^(0)]  = (0.75, 3/2;-0. 350,0)  and  tR/t1  = I.06  with 
a -5$  second-order  drift  correction;  the  center-vortex  trajectory 
(and  envelope)  is  reminiscent  of  a large  center-vortex  initial 
displacement  in  the  hard-vortex  regime  (cf.  Fig.  15);  in  fact  the 
motion  of  the  center  vortex  is  outside  of  the  center  stability 
boundary  shown  in  Fig.  18,  albeit  for  relatively  short  periods  of 
time;  this  behavior  differs  from  the  hard-vortex  regime  where  the 
center-vortex  trajectory  is  always  confined  by  the  stability 
boundary.  In  Fig.  22,  [>,n;  x^1  ^ ( 0), y^ 1 ^ ( 0) ] = (0.5, 5/4; 0. 264, 0) 

A 

and  tR/t-^  = 1.25  with  a -1.5%  second-order  drift  correction;  for 
positive  initial  radial  displacements  the  vortex  trajectories  are 
confined  within  their  respective  stability  boundaries;  and  each 
vortex  makes  four  excursions  around  its  equilibrium  position.  Two 
neighboring  recurrent  trajectories  are  shown  for  k = 0.25:  In 

Fig.  23a,  [fc,n;x^1^(0),y^1^(0)]  = ( 0. 25,  25/19;  0.452,  0)  and 
tR/ti  = I.58  with  a -3.2%  second-order  drift  correction  and  19 
excursions  in  a recurrence  time;  in  Fig.  23b,  [jc^x^1  ^ ( 0), yj x ^ ( 0) ] 

= (0.25,4/3; 0.456, 0)  and  tR/t^  = 1.6l  with  a -3. 3%  second-order 


drift  correction  and  only  three  excursions  in  a recurrence  time; 
the  difference  in  the  initial  displacement  is  just  0.004  and  the 
difference  in  the  frequency  ratio  An  = 1/57  between  the  two  runs; 
this  is  the  best  resolution  we  have  made  of  neighboring  frequency 
ratios.  In  Fig.  24a,  U,n;x^(0),y^  ( 0) ] = ( 0. 125,  27/26; 0.250,  0) 

A 

and  tR/t^  = 1.20  with  a -1.15$  second-order  drift  correction  and  26 
excursions  in  one  recurrence  time;  without  this  small  second-order 
correction,  a considerable  drifting  motion  occurs  as  shown  in 
Fig.  24b.  Three  recurrent  trajectories  are  shown  for  k = 0:  In 

Fig.  25a,b,c,d  (/c,n)  = (0,1)  with  0.001  initial  perturbation  for 

a)  ^"^(O),  b)  y^^O),  c)  x^1)  ( 0)  and  d)  y^(0);  negative- 
perturbation  trajectories  are  obtained  by  rotating  the  positive- 
perturbation  trajectories  by  tt  around  the  equilibrium  position  for 
each  of  the  vortices;  the  initially  displaced  center-vortex  motion 
is  a linear  oscillation  (Fig.  25c, d);  in  Fig.  26a, 
[«:,n;x^1)(0),y^l)(0)]  = ( 0, 8/7;  0.503, 0)  and  tR/t±  = 1.53  with  a 
-3.53$  second-order  drift  correction  and  seven  excursions  in  one 
recurrence  time;  in  Fig.  26b,  [/c, n; x^1^  ( 0), y^1^  ( 0)] 

= (0, 6/5; 0.300, 0)  and  tR/t^  = 1.09  with  a -2$  second-order  drift 
correction;  in  all  these  cases  the  trajectories  are  confined 
within  their  stability  boundaries.  Finally  two  unstable  trajec- 
tories are  shown  for  k - 0:  In  Fig.  27a,  [x^'*‘^(0),y^1^  (0)] 

= (-0.300,0)  with  7q  = -1.0  and  no  drift  correction  whatsoever  in 
which  case  the  center- vortex  is  absolutely  unstable  and  goes  to 
infinity4 with  one  of  the  circle  vortices  while  the  remaining  two 
vortices  are  left  to  themselves  in  a circular  trajectory;  in 
Fig.  27b,  [x^(0),y^(0)]  = (0.600,0)  and  = -0.93500  with  a 
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-6.5#  second-order  drift  correction;  in  contrast  to  the  run  with- 
out a drift  correction,  the  center- vortex  is  soon  captured  by  one 
of  the  circle  vortices  for  a nearly  circular  excursion  far  outside 
the  original  circle,  then  goes  on  to  a second  trip  with  the  next 
circle  vortex,  progressing  in  a counter-clockwise  order. 

2.  The  Instability  or  Critical  Point,  >ccr  = 1.44763516  (n=2) 

The  problem  of  determining  by  numerical  calculation  that 
this  (3+l)-vortex  configuration  is  unstable  at  kct  is  not  easy 
since  the  configuration  apparently  tends  toward  stability  in  an 
algebraic  manner  as  the  initial  perturbation  is  decreased  toward 
zero;  in  this  particular  circumstance,  a second-order  analytical 
calculation  is  being  carried  out,  which  may  turn  out  to  be  easier 
and  will  verify  our  present  results.  In  the  following  sequence  of 
runs  we  find  that  certain  types  of  equal-amplitude  initial  dis- 
placements are  more  stable  than  others;  this  structural  feature 
compounds  the  numerical  difficulties.  In  Fig.  28a,b,c,d, e the 
initial  displacements  are  a)  x^^(0)  = -0.005#  b ) y^^(O)  = 0.005, 
c)  y^(0)  = -0.005,  d)  x^(0)  = 0.005  and  e)  yj^(0)  = 0.005; 
all  other  component  displacements  are  in  their  equilibrium  posi- 
tions in  all  of  these  figures.  Only  in  Fig.  28e  is  instability  an 
uncertainty  and  those  initial  displacements  not  shown,  viz. 
x^^(0)  = 0.005  and  x^^(0)  = -0.005#  are  just  as  unstable  as 
those  of  opposite  sign  (cf.  Figs.  28a  and  2dd).  Thus  in  Fig. 

29a, b are  shown  the  trajectories  for  a)  y^'(0)  = 0.015  and 
b)  y^(0)  = 0.02;  the  center-vortex  trajectory  in  Fig.  29b  is 
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clearly  becoming  unstable.  Finally  in  Fig.  30a,b  for  a) 
[^(Ohy^^O)]  = (0.01,0.03)  and  b)  LX^1)(0),y^)1)  (0)] 

= (-0.01,0.03),  the  center  vortex  goes  outside  of  the  circle  and 
is  captured  by  one  or  more  circle  vortices  in  sequence. 

3 • Unstable  Motions 

Our  study  of  unstable  motions  is  limited  in  scope  because  of 
our  self-imposed  restriction  to  emphasize  those  motions  in  which 
the  center  (cyclonic)  vortex  remains  inside  the  circle  of  anti- 
cyclonic  vortices  and  where  the  vortex-strength  ratio  is  such 
that  the  stationarity  condition  holds.  These  motions  are  unstable 
only  in  the  sense  that  the  trajectories  do  not  remain  in  the 
immediate  neighborhood  of  the  stationary  equilibrium  positions  for 
each  of  the  four  vortices;  this  is  a necessary  condition  from  the 
viewpoint  of  the  climatological  problem  but  possibly  not  for  other 
physical  applications  of  these  hard  and  soft  vortices.  Otherwise 
we  could  regard  these  motions  as  not  unstable  but  as  a change  in 
mode  from  motions  near  one  equilibrium  state  into  some  other 
quasi-equilibrium  motions.  Here  we  show  two  more  examples  of 
unstable  motions  and  a summarizing  result  concerning  the  strength 
of  the  nonlinear  interaction  between  vortex  pairs  moving  in  the 
second  mode  of  quasi-periodic  motions;  these  results  include  many 
runs  (more  than  600  or  approximately  10i&  of  the  total  number  of 
runs)  over  the  entire  range  of  the  length  parameter  0 < k < 5.  In 
Fig.  31  unstable  trajectories  are  shown  for  k = 5 and 
[Xq1  ^ (0),y£ 1 ^(0)]  = (0.210,0);  the  order  of  capture  of  the  center 
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vortex  by  one  of  the  circle  vortices  is  2-3-1-3-2-1-3  which  shows 
that,  after  the  first  two  irregular  transfers  the  capture-ordering 
appears  to  progress  in  a regular  (hard-vortex)  clockwise  direc- 
tion; and  the  center  and  circle  vortices  progress  in  a clockwise 
direction  with  nearly  circular  loops  in  phase  with  each  other 
during  the  time  of  capture.  In  Fig.  32  unstable  trajectories  are 
shown  for  k = 2.5  and  [x^1^  (Ohy^1^  (0)]  = (-0.130,0);  here  the 
circle  vortex  is  captured  for  three  transits  by  the  #1  circle 
vortex. 

In  the  case  of  climatological  instability  the  vortex, 
originally  at  the  center  of  the  circle,  leaves  the  vicinity  and 
follows  a path  going  around  first  one  of  the  circle  vertices,  then 
another  one,  etc.  Because  of  the  nonlinearity  of  the  problem  it 
is  virtually  impossible  to  predict  in  what  order  the  "center" 
vortex  will  be  captured  by  the  circle  vortices.  Certain  general- 
ized observations  about  these  trajectories  can  nevertheless  be 
made:  When  the  motion  is  still  in  or  near  the  region  of  stability, 

the  center  vortex  will  spend  some  time  — several  loops  — near  the 
center.  Then  it  takes  off  suddenly  (exponentially),  encircles  one 
of  the  circle  vortices,  and  returns  to  the  center,  only  to  take 
off  again  later  on.  In  this  case  the  order  in  which  the  circle 
vortices  are  taken  up  is  completely  unpredictable  numerically. 
Farther  out  in  the  region  of  instability,  the  center  vortex  takes 
off  almost  immediately,  and  after  circling  one  of  the  circle 
vortices,  takes  up  another  one,  etc.,  often  in  a systematic  order. 
Very  occasionally  it  will  circle  one  of  the  outer  vortices  twice 
or  mere  times  in  succession.  For  soft  vortices  there  also  exists 
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the  possibility  of  total  capture  by  only  one  circle  vortex,  i.e. 
the  center  vortex  encircles  one  of  the  circle  vortices  almost 
indefinitely,  while  the  remaining  two  vortices  circle  each  other 


The  radii  of  the  approximate  circles  formed  by  the  trajec- 
tories appear  to  be  of  some  significance,  and  bear  some  relation- 
ship to  the  value  of  the  vortex  strength  ratio  yo.  In  order  to 
make  a good  numerical  approximation  to  the  radii  of  these  circles, 

the  numerical  results  were  analyzed  as  follows:  Whenever  the 

2 2 

coordinates  of  the  ''center"  vortex  become  such  that  x +y  >1, 
i.e.  whenever  that  vortex  left  the  interior  of  the  unit  circle, 
the  coordinates  at  all  succeeding  time  intervals  were  stored  in 
the  computor  memory  until  the  vortex  returned  to  the  interior  of 
the  circle.  The  coordinates  so  obtained  were  fitted  by  a least- 
square  procedure  into  a circle,  and  the  center  and  radius  of  that 
circle  computed.  The  same  calculation  was  done  at  the  same  time 
for  the  coordinates  of  the  circle  vortices.  An  inspection  of  the 
results  for  the  four  vortices  made  it  a simple  matter  to  decide 
which  two  had  a common  center  (one  vortex  going  in  phase  around 
the  other).  The  radii  of  these  two  vortices  were  then  determined, 
and  the  ratio  of  these  radii  were  found  to  be  closely  correlated 
to  I7  In  fact  in  the  absence  of  the  other  two  vortices,  this 

Radii  Ratio  is  equal  to  the  inverse  ratio  of  the  vortex  strengths. 
In  Fig.  33  the  variation  of  | -y Q f • (Radii  Ratio)  with  respect  to  /c 
is  shown  for  0 < ic.  ^ 5:  1)  it  is  nominally  equal  to  one  for  k = 0; 

2)  it  decreases  to  a minimum  value  equal  to  0.7  at  < = 4;  and 

3)  it  apparently  increases  for  k.  > 4 (this  increase  may  be  attrib- 
utable to  numerical  inaccuracies  for  (<r)  > 10). 


In  Ref.  35,  Appendix  B the  sum  E = ^ ^krk0k  has  ^een  s^own 
to  be  constant  for  1)  k.  = 0 and  2)  after  linearization;  we  have 
pointed  out  in  Footnote  5 that  E is  usually  not  an  invariant  for 
k > 0.  This  sum  may  be  useful  for  studying  the  unstable  motions, 
as  well  as  for  future  studies  cf  the  statistical  properties  of 
Bessel  vortices.  Compared  to  logarithmic  vortices  (k  = 0),  these 
Bessel  vortices  have  better  physical  properties:  1)  finite 

kinetic  energy  and  2)  one  less  integral  constraint. 
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V.  Remarks 


A.  An  Interpretation  of  Soft  and  Hard  Vortex  Stability- 

In  Ref.  3 the  applicability  of  the  study  of  kinematic  (quasi- 
dynamically-balanced,  time-asymptotic)  motion  to  a large  range  of 
length  scales  in  suitable  physical  problems  has  been  pointed  out. 
Our  present  study  supports  this  viewpoint  for  at  least  two  differ- 
ent length  (and  time)  scales  of  atmospheric  motions:  1)  In  Ref.  9 

our  barotropic  model  has  been  applied  to  the  tracking  of  hurricane 
motion,  using  one  Bessel  vortex,  where  the  normalized  length  par- 
ameter k was  chosen  in  the  range  1 _<  < < -4/3;  this  <-range  in  our 
present  context  is  in  the  soft-vortex  regime;  while  2)  our  (3+1)- 
vortex  representation  of  global-scale  hemispherical  motions  is  in 
the  hard-vortex  regime.  One  possible  interpretation  of  our  stabil- 
ity results,  showing  the  existence  of  two  distinct  stable  regimes 
separated  by  a critical  instability  point  kcy  (with  n = = 2), 

is  that  there  is  a short-range  moderate-scale  atmospheric  descrip- 
tion available  as  well  as  a longer-range  global-scale,  kinematic 
representation;  but  an  intermediate-scale  regime  of  relatively  dis- 
ordered motions  exists  where  some  statistical  information  may  be 
necessary.  This  conforms  to  our  qualitative  observations  of 
natural  phenomena  in  which  there  is  a sequence  of  ordered  and  dis- 
ordered motions  depending  on  different  length  and  time  scales.  For 
the  short-range  description  of  atmospheric  motions,  there  is  appar- 
ently an  upper  limit  of  about  three  or  four  weeks  to  the  predicta- 
bility of  the  physical  variables  using  a kinematical  (or  dynamical) 
model;  and  the  use  of  combined  kinematical-statistical  model  has 
been  studied  recently  by  a number  of  atmospheric  scientists,  e.g. 
Leith  [10]. 
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B.  Applicability  of  the  (3+l)-Vortex  Configuration  to  Climatology 
It  may  seem  surprising  that  this  simple  configuration  of 


four  K^-vortices  can  be  at  all  representative  of  atmospheric 
motions.  Nevertheless  our  present  study  shows  the  apparent  valid- 
ity of  this  planar  model  in  simulating  some  basic  properties  of 
global-scale  motions  with  correspondingly  long  time-scales  greater 
than  a season,  a year  or  tens  of  years.  One  factor  which  gives  us 
some  confidence  in  the  reality  of  this  lowest  order  barotropic 
model  is  that,  in  the  climatological  range  of  the  length  parameter, 
the  vortex  strength  ratio  is  small  (less  than  0.2). 

A few  runs  were  made  to  estimate  the  effect  of  the  variation 
of  the  Coriolis  parfimeter  between  the  anticyclonic  latitude  circle 
and  the  pole.  Both  k (circle  vortex)  and  kq  (center  vortex)  were 
varied  in  the  range  2 <•  (k;/cq)  < 4.  We  find  that  the  value  of  ic 
essentially  controls  the  stability  properties  over  a wide  range  of 
values  of  tcQ.  For  example  comparing  (k;k0)  = (4;2)  with  k=kq=^, 
the  stability  boundaries  for  circle-  or  center-vortex  initial  dis- 
placements differ  by  less  than  0.01. 

The  extension  of  this  approach  toward  more  realistic  geometry 
and  configurations  is  being  pursued.  Peters  [11],  and  Leiva  [12] 
have  recently  begun  the  study  of  vortex  motions  on  a rotating 
sphere.  Leiva* s dissertation  shows  some  of  the  quantitative 
connections  between  the  linearized  stability  analysis  on  the  sphere 
as  compared  to  the  planar  approximation;  the  qualitative  relation- 
ships are  quite  similar;  however  a more  complete  study  of  the  non- 
linear motions  is  required. 
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Neutral  Stability  Diagn 


70  vs.  K for  O = 0 and  N = 2,  3 and  4 
Stability  Curves  for  N = 2,  3 and  4. 


and  Neutral 


Initial  Displacement 


Fig.  4 Stability  Boundaries  for  (N=3;O=0)  and  Initial  Displacement 
x^^(0)  of  the  Center  Vortex  for  0 < k < 5. 


al  Displacements  [jq  ;(0),yi  1M0)1  of 
m m i,  j 3 

x'  1 ,y\  '(0)]  of  the  Center  Vortex 


Pig.  6 Range  of  Stable  Initial  Displacements  1 > (0),y< (U (0 j ] 

of  a Circle  Vortex  for  a:  =3.2. 


Fig.  8 Range  of  Stable  Initial  Displacements  (0),y^^(0)] 

of  the  Center  Vortex  for  k = 4. 
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Vortex  Trajectories  for  « = ? and  Small  Initial  Displacements: 
(a)  x^1}(0)  = 0.001,  (b)  x^lj(0)  = -0.001,  (c)  y^fo)  = 0.001 
and  (dj  x^^(0)  = 0.001;  Magnification  500X. 
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Fig.  11  Vortex  Strength  Ratio  I^J"1  vs.  n > 2 l'or  < = 1.75,  2.5, 
3,  3»5>  3. 75*  ^ and  5. 


Vortex  Trajectories  for 
fa)  [<c, njx^1  ^ (0),y^1  ^ ( 0 

( b ) k = 3.2  with  [ X4 ^ ( 1 
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(<i)  k = 3.2  with  [x. 
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k - 3.2  with  [x. 


The  second  order  drift  corrections  are 
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Fig.  18  Range  of  Stable  Initial  Displacements  [x^^(0),y^^(0)]  of  a Circle 
Vortex  and  [x  1^(0),y^1^(0)]  of  the  Center  Vortex  for  *■  = 0.75  . 
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(a)  Vortex  Trajectories  for  [<c,n 

= (0. 25, 25/19;0. 452,0). 


(b)  Vortex  Trajectories  for  [«,n;x. 

= (0. 25, 4/5;0. 456,0). 


Vortex  Trajectories  for  (/c,n)  = (0, 1)  with  0.001  Initial 
Displacements:  (a)  xJ^'fO),  (b)  vi1^(0),  (c)  x'^(0)  and 

(d)  y'  ;(0);  Magnification  500X. 


Fig  27o 


Fig  27b 

(a)  Vortex  Trajectories  for  <e  = 0 and  [x^ 1 ^ (0 ), y 1 ^ ( 0)  ] 

= (-0.300,0);  7o  = -1;  no  drift  correction;  absolutely 

unstable . 

(b)  Vortex  Trajectories  for  <c  - 0 and  [x^  1 ^ ( 0 ), yi  1 ^ (0 ) ] 

= (0.600, 0);  y,  = -0.93^00;  unstable. 
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Vortex  Trajectories  for  k.  = 1.4476  and  Initial 


Displacements : 
y^1}(0)  = -0.005, 


Unstable  Vortex  Trajectories  for  «c  = 2.5  ana 
[x^1)(0),y^1,(0)]  = (-0-130,0). 
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